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Abstract 

Published measurements of semi-inclusive A production in p-Au collisions at the 
AGS are used to estimate the yields of singly strange hadrons in nucleus-nucleus 
(A-A) collisions. Results of a described extrapolation technique are shown and com- 
pared to measurements of K + production in Si-Al, Si-Au, and Au-Au collisions at 
the AGS and net A production in S-S, S-Ag, Pb-Pb, and inclusive p-A collisions at 
the SPS. The extrapolations can account for more than 75% of the measured strange 
particle yields in all of the studied systems except for very central Au+Au collisions 
at the AGS where RQMD comparisons suggest large re-scattering contributions. 



Strange particle production has long been considered an important experimen- 
tal probe of heavy ion collisions due to the possibility that strange particle 



Preprint submitted to Elsevier Preprint 



8 February 2008 



yields may be significantly enhanced by quark-gluon plasma (QGP) forma- 
tion [1,2]. Heavy ion experiments at the Brookhaven National Laboratory 
AGS and CERN SPS accelerators have observed factors of 3-4 enhancements 
in the yield per participant of strange mesons and singly strange baryons [3- 
10] and, at the SPS, larger enhancements in the production of multiply strange 
baryons [11-13]. Cascade models can quantitatively reproduce the enhance- 
ment of strange mesons in central heavy ion collisions at the AGS but are 
unable to reproduce the enhancements observed at the SPS without the intro- 
duction of "exotic" processes. Furthermore, inclusive measurements of kaon 
production in proton-nucleus (p-A) collisions at the AGS have shown an en- 
hancement with target size [14,15], while no such enhancement was seen at 
the SPS [16]. The strong enhancements in nucleus-nucleus (A-A) collisions at 
the SPS, combined with the lack of enhancement in p-A and the failure of cas- 
cade models, have been cited as experimental evidence for QGP formation at 
the SPS [17,18], while the enhancements at the AGS are typically attributed 
to rescattering of produced hadrons. Then, the similarity in magnitude and 
pattern of enhancements of singly strange hadrons at the AGS and SPS must 
be attributed to coincidence. 

However, previous measurements of semi-inclusive A production by Experi- 
ment 910 at the AGS [19] suggest an alternative, common explanation for 
the enhancement of singly strange hadron production in A-A collisions at the 
AGS and SPS. E910 has measured the dependence of the total A yield in 
17.5 GeV/c p-Au collisions on the number of scatterings v of the proton in 
the Au nucleus and demonstrated enhanced A production over a number of 
participant (Apart) scaling of p-p data. In this paper we present the results of 
an analysis that extrapolates the effect observed by E910 to A-A collisions and 
predicts the hyperon-associated K + yields in Si-Al, Si-Au, and Au-Au colli- 
sions at the AGS and the net A yields in S-S, S-Ag, and Pb-Pb collisions at the 
SPS. We show that the extrapolation of the effect observed in p-A collisions 
can account for most of the observed enhancements in hyperon-associated K + 
production at the AGS and A production at SPS in both light and heavy ion 
induced reactions. With the same analysis we also explain the lack of observed 
enhancement in inclusive p-A measurements at the SPS. 

E910 reported [19] that the A excess in p-Au collisions over a simple A part or 
wounded-nucleon scaling of p-p data increased proportionally to v for the first 
three collisions with a slope equal to one half of the p-p A multiplicity, 

AJV A (i/) = (i/-l)iJVf. (1) 



Since the projectile contributes |A^ P to the total A multiplicity in the first 
collision, the E910 excess is consistent with a projectile contribution to the 



total A multiplicity that grows with v, as [19] 

AT° j = \N™ u,u<3, 



(2) 



and remains constant for v > 3. To extrapolate to A-A collisions we assume 
that Eq. 2 is valid for all relevant beam energies and that it applies equally to 
all participating nucleons with A^ p replaced by the estimated yield in isospin 
averaged nucleon-nucleon (N-N) collisions, N NN . Then, the total multiplicity 
in an A-B collision is 

N AB = 1 N NN £ [A p M + B p t{v)] y + | N NN £ APp{u) + B p t{u) (3) 

P p /t(^) represents the probability distribution for the projectile and target nu- 
cleons to undergo v primary scatterings during the collision. Since the number 
of participants can also be calculated from the P p / t (z/), 

N pi , It AB = J2AP p (v)+BP t (is), (4) 

we can directly obtain the relationship between N AB and N paTt - We note that 
Eq. 3 applied to p-A collisions reduces to 

A^pA = N NN z/, v < 3 ; N pA = -N NN (z/ + 3), v > 3. (5) 

The total A multiplicities measured by E910 saturate for v > 5 implying 
a violation of Eq. 5. This saturation was attributed to the stopping of the 
projectile [19] and is not expected to be present at SPS energies. 

In our analysis the P p / t (z/) distributions were calculated using both a Glauber 
model [20] and the Lund Monte-Carlo [21]. The results are presented and 
compared in Fig. 1 in terms of R part = 2N AB /(A part N NN ), i.e. the factor by 
which the A-A yield per participant is enhanced relative to N-N collisions. In 
the Glauber calculation, the P p / t (^) distributions were evaluated at fixed val- 
ues of collision impact parameter b, and the simultaneous dependence of N AB 
and A" par t on b was used to obtain N AB (Ap art ). This procedure, while calcu- 
lationally simple, potentially introduces errors by averaging over fluctuations 
that might differently affect N AB and Ap art . The slight differences between 
the two calculations at large b are consistent with the expectation that the 
Monte-Carlo calculation would better account for these fluctuations. The two 
calculations give very similar results for the A" part dependence of R part , though 
the Monte-Carlo result has a slightly slower growth in R part with A" part . To 



evaluate the importance of the stopping effects at AGS energies implied by 
the E910 data, we additionally imposed in the Glauber calculations a max- 
imum thickness to the nuclei corresponding to 5 nucleon interaction lengths 
with the result shown in Fig. f by the dashed lines. Hereafter we will use the 
Monte-Carlo results when presenting the calculation results. 

We now compare our calculations to experimental data starting with K + mul- 
tiplicities measured by Experiments 859 and 866 in Si-Al, Si-Au, and Au-Au 
collisions at the AGS [4]. At AGS energies, the K + is predominantly produced 
in association with hyperons, so the total K + multiplicity should have the 
same v dependence as the A (excluding K + K contributions). Previously esti- 
mated production rates for K + production in isospin-averaged N-N collisions 
are listed in Table 1 for the three different colliding systems [4] . At comparable 
energies p-p measurements [22] suggest that 15% of the K + yield results from 
K + K processes. Lacking clear knowledge of how these processes grow in p-A 
collisions we conservatively assumed that K + K production increases propor- 
tionally to iVp art . Thus, the N-N K + rates used in Eq. 2 were reduced by 15%, 
and the K + K contribution to the K + yield was added separately. Figure 2 
shows a comparison between the results of our calculation and the experimen- 
tally measured K + yields. Two sets of points are shown for the Si-Au data 
because the experimentally quoted A part values were found to be inconsistent 
with calculated values of A part for the quoted centrality fractions [4], though 
nowhere by more than twice the stated error on the A part values. Since the two 
calculations give results that differ by less than 2%, we attribute the discrep- 
ancy with the data to a systematic error in the method used by the experiment 
to estimate A part in asymmetric collisions. The shaded region in Fig. 2 indi- 
cates the uncertainties in our extrapolation due to the systematic errors on 
the N-N yields shown in Table 1 and the uncertainty in the K + K rate in p-p 
collisions. Our calculations account for more than 75% of the measured K + 
yields everywhere except in the most central Au-Au collisions where the data 
exceeds our extrapolation by 50%. We also show in Fig. 2 the results from 
the cascade model RQMD [23] simulating 11.1 AGeV/c Au+Au collisions. 
Although RQMD significantly under-predicts the E866 K + yield over most of 
the centrality range, it shows a rapid increase in K + yield with iV part for very 
central collisions that may be due to rescattering processes, which would be 
expected to grow roughly like A part 2 . The rapid rise in the RQMD calculation 
is qualitatively similar in shape to the observed excess of the central Au-Au 
data over our calculation suggesting that this excess may, in fact, be due to 
rescattering processes not present in p-A data. 

We also apply our extrapolation procedure to SPS energies to determine 
whether the strangeness enhancements there may also be accounted for by the 
effect observed in the E910 data. Since Eq. 2 may not describe baryon-anti- 
baron (BB) production which is negligible at AGS energies but important at 
SPS energies, comparisons with SPS data are done using the "net" (N A — N^) 



yields. Figure 3 shows measurements of total net A production from NA35 
in central S-S [8] and S-Ag [10] collisions at 200 A-GeV/c and net A produc- 
tion at mid-rapidity in Pb-Pb collisions at 160 A-GeV/c from WA97 [12]. The 
input N-N rates required for the extrapolation are listed in Table 1. For the 
Pb-Pb extrapolation we reduced the p-Be A yields of WA97 by a factor of 0.87 
to remove the enhancement that would result from multiple collisions of the 
proton according to our method. The extrapolations, shown in Fig. 3 by solid 
lines, reproduce the observed A enhancement in S-A collisions but fall well be- 
low the Pb-Pb data. However, the problem with the latter comparison is that 
the WA97 measurements are limited to an acceptance of unit rapidity about 
the N-N center of mass. The p-p A rapidity distribution is strongly peaked 
at forward/backward rapidity and suppressed near mid-rapidity [24,25]. The 
increased stopping of the baryons in Pb-Pb collisions [26] will strongly shift 
the peaks toward mid-rapidity, so that an enhancement would be observed 
by WA97 even if there were no enhancement in the total A yield. From pa- 
rameterizations of the A xp distribution in p-Be collisions [27] we estimate 
that 11% of the total p-Be A yield is contained within the central unit of 
rapidity. Using recent measurements and extrapolations for A and A dn/dy 
for 158 A-GeV/c Au-Au [28] we estimate that 21% of the total net A yield is 
contained in the central unit of rapidity. If we multiply our extrapolation by a 
factor of 0.21/0.11 = 1.9 to account for the effects of the dn/dy shape change, 
we then obtain the dashed line shown in Fig. 3. The excellent agreement be- 
tween this result and the Pb-Pb data may be partially accidental since we 
have not included an N paTt dependence to the correction factor. Nonetheless, 
we can conclude that our extrapolation accounts for more than 75% of the net 
hyperon yields in light and heavy ion collisions at the SPS. 

The large enhancements at SPS energies predicted by the E910 extrapolation 
appear to contradict the lack of enhancement noted by previous p-A mea- 
surements at the SPS [16]. However, the effect observed by E910 may be sub- 
stantially obscured in inclusive collisions due to the large contributions from 
peripheral collisions and due to the saturation of the projectile enhancement 
at v — 3, which will result in a decreasing A yield per participant for v > 3. 
Figure 4 shows results of our calculations for inclusive collisions compared 
with a variety of inclusive measurements [16]. Our calculation reproduces the 
nearly constant A yield per participant observed in the experimental data in 
spite of the enhanced projectile contribution. 

In conclusion, we have performed a set of calculations intended to quantita- 
tively evaluate the significance of the observed enhancement of A production 
in p-Au collisions [19] in the interpretation of strangeness enhancement in 
heavy ion collisions. Our extrapolation is based on applying the enhancement 
observed by E910 to all participant nucleons in A-A collisions assuming that it 
results directly from the multiple interactions of the proton. We show that our 
extrapolation can account for more than 75% of the total hyperon or hyperon 



associated kaon yield in light and heavy ion induced collisions at the AGS and 
SPS except for K + production in central Au-Au collisions. Here, comparisons 
with RQMD suggests that the discrepancy is due to rescattering processes 
that contribute most importantly in central collisions. The analysis presented 
here suggests that a dynamical mechanism and not QGP formation may be 
responsible for enhancing the yields of (at least) singly strange hadrons. Our 
analysis also suggests that the similar enhancements observed at the AGS 
and SPS result from a common mechanism in spite of the difference in beam 
energies. 
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Fig. 1. Calculated enhancement factor R par t (see text for details) vs b and A part for 
various colliding systems, solid lines - Glauber geometry, points - Lund geometry. 

Fig. 2. Comparison between measured K + yields in Si-A (left) and Au-Au colli- 
sions (right) at the AGS and extrapolations, ■ - Si-Al, A - Si-Au, • - Au-Au, T - 
Si-Au with adjusted Apart, solid line - extrapolations, dashed - A part scaling of N-N, 
dot-dashed - RQMD calculation for Au-Au. Shaded regions indicate la systematic 
errors due to uncertainties in N-N yields. 



Fig. 3. Comparison of net A (A — A) yields from (left) S-A and (right) Pb-Pb col- 
lisions at the SPS and calculations, T - S-S, ▲ - S-Ag, • - Pb-Pb. Lines, solid 
- extrapolations, dashed - A^ part scaling of N-N, dot-dashed (right) extrapolation 
corrected for stopping (see text for details) 

Fig. 4. Comparison of average net A yield per participant in inclusive p-A colli- 
sions 200 GeV/c and extrapolation of E910 data, • - p-A, □ - p-p, solid line - 
extrapolation. 
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